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^ , The electrostatic perturbations in an unmagnetized, non-isothermal [Ti<^Te) 

H : 

electron-ion plasma shear flow are considered. New physical eff'ects, arising 
due to the non-normality of linear dynamics, are described. It is shown that 
the velocity shear induces the extraction of the mean flow energy by the acous- 
tic perturbations (ion-sound waves). The influence of the medium dispersion, 
rising due to the violation of quasineutrality for perturbations, is examined. 
It is shown that in the course of the evolution ion-sound waves turn into 
ion plasma oscillations. New class of nonperiodic, electrostatic perturbations 
(with vortical motion of ion component), characterized by the intense energy 
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exchange with the mean flow, is also described. 
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Classical stability theory of continuous media motion (normal mode approach) has been 
victorious in explaining how different kinds of shear flows become unstable. However, in 
some, quite simple and important kinds of parallel free shear flows (e.g., plane Couette and 
Poiseuille, or pipe Poiseuille flows) the approach has serious problems, evoked by the not 
self-adjoint character of the governing equations [1-3]. That is why the predictions of the 
traditional stability approach fail to match the results of most experiments with these kinds 
of flows [3] . 

An alternative approach to the problem is that due to Kelvin [4], which implicates the 
change of independent variables from a laboratory to a moving frame and study of the non- 
exponential temporal evolution of spatial Fourier harmonics (SFH). The method is operative 
for any mean velocity proflles but it is most manageable to ones that are linear, or peacewise 
linear [1,5]. Effectiveness of the method was repeatedly proved, when it helped to obtain 
the bunch of unlooked-for results on the dynamics of the perturbations in hydrodynamical 
[6-12] and hydromagnetic [13-15] shear flows. 

In [12] authors have considered the evolution of two-dimensional (2D) SFH in a com- 
pressible, plane hydrodynamical Couette flow. Their analysis, which involved the nonmodal 
approach, brought them to the discovery of the new mechanism of the energy exchange be- 
tween the mean flow and sound-type perturbations. In particular, it was shown that energy 
of SFH may grow linearly in time — perturbations extract the energy from the mean shear 
flow. This process appears to be quite universal and one should expect that it may be also 
influential in the wide variety of continuous media with the analogous kinematics. 

In this report we shall examine the case of electron- ion plasma shear flow and show that 
the process of velocity shear induced energy transfer from the mean flow to the collective 
modes exists and can be quite efficient in this case too. Moreover, as we shall see, the 
peculiarity of the plasma state of the medium plays its distinctive role and leads to the 
whole group of interesting new effects. 

As it is well-known, in the coUisionless unmagnetized plasma with Tj<^Te (where Tg and 
Tj are electron and ion temperatures respectively) exists the weakly damped low-frequency, 
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longitudinal, and therefore electrostatic, ion mode. When its wave-length A~l//c {k is 
wave number) greatly exceeds electron Debye length (A^A£)e=(T'e/47re^no)^/^, then this low- 
frequency mode represents the ion-sound wave with constant phase velocity ('ion acoustic' 
speed): Cs={Te/miY^'^ . However, the velocity shear induces the "linear drift" of SFH (the 
process well- acknowledged in the literature [6,12,13] cultivating the nonmodal approach) 
which, mathematically, is exposed in the temporal variation of the wave number vector k(t). 
It means that the influence of the dispersion of ion mode, arising as the result of the violation 
of quasineutrality for the perturbations, which, in other words, is related to the finiteness of 
kXoe-i should be taken into due account. The reason is simple: if initially \k\De\^^ the linear 
drift process will, eventually, transfer the SFH to the region of k — space, where the latter 
condition does not hold. In physical terms it means that that under certain circumstances 
the ion sound waves, drawing energy from the mean shear flow, subsequently turn into the ion 
plasma oscillations. The latter collective mode is weakly damped if A»ADi=(Ti/47re^no)^/^. 

Velocity shear of the mean flow induces, also, excitation of the completely new kind of 
non-periodic, electrostatic perturbations with vortical motion of the plasma ion component. 
These perturbations are able to the effective exchange of their energy with the mean flow 
and under certain conditions may play the dominant role in the behavior of the plasma flow. 

At the end of this report we shall point out on the possible relevance of above noted re- 
markable processes to the problems of excitation of electrostatic waves, resonant acceleration 
of particles and onset of a turbulence in non magnetized electron-ion plasmas. 

Let us consider electrostatic oscillation modes in the shear flow of electron-ion nonmag- 
netized plasma. The mean velocity Uo=(>lt/, 0, 0) is directed along the X axis and has a 
hnear shear along the Y axis. Without the loss of generahty we can assume that the con- 
stant parameter A is positive. Equilibrium electric field E is zero. Taking approval of the 
electrostatic nature of oscillations we can write the electric field perturbation as the gradient 
of the electric potential: E = — V0. The plasma in its equilibrium state is assumed to be 
homogeneous and quasineutral — equilibrium number density of electrons n^Q is equal to the 
one of ions uiq and is constant: n^o — rijo — uq — const. 
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When the plasma is properly rarefied, ion and electron components interact weakly with 
each other and may, thus, have very different temperatures. In particular, we shall assume 
that temperature of electrons is homogeneous Tg = const and for the sake of simplicity take 
that Ti — 0. The effects, which may be evoked by the finiteness of Ti will be discussed briefly 
at the end of the report. 

The usual approximation for electrostatic low-frequency waves implies that the electron 
number density can be described by the Boltzmann distribution [16]: 

ne = noexp(e0/re)Rino(l + e^/Te), (1) 

Let us, now, decompose the instantaneous values of all physical variables onto their mean 
and perturbed components: Vj = Uq + u, ne,i — uq + n'^ .^. The basic system of linearized 
equations for ions, describing the evolution of small-scale, 3D perturbations in this flow, 
takes the form: 

[dt + Ayd^jn'i + no[da:U^ + dyUy + d^^u^] = 0, (2) 
[dt + Ayd,]u + (u, V)Uo = -— V0, (3) 

m 

'9' + 9.' + 9']0 = 47reK-nO. (4) 

where (u, V)Uo is a vector, which has only one, x-th, nonzero component equal to Auy. 

Let us make the following substitution of variables: Xi — x—Ayt; yi — y; Z\ — z\ t\—t 
and perform the Fourier analysis of (2-4), expanding unknown functions with respect to 
spatial coordinates: 

F = y" dkxidky^dkz^F{kxi, ky^, k^^, ti)exp[i{kxiXi + ky^yi + k^^zi)], (5) 

where under F we imply all physical variables appearing in the problem. Instead of equations 
(2) and (3) we get (for details, see [9,13]) the following set of first order, ordinary differential 
equations (ODE's) for SFH of these variables: 

7Vf)=^;, + /3(T)^, + 7^., (6) 
^(^) = -Rvy - (7a) 



= -/3(t)$, (76) 
= -7$; (7c) 
complemented by the algebraic relations, arising from (1) and (4): 

N, = (8) 

iV,= [l + e'/C^JiVe, (9) 

where, hereafter, F^") will denote the n-th order time derivative of F. 

In (6-9) we have introduced the following dimensionless notations: R=A/Cskxi, 
T=Cskxiti, PQ=kyj^/kxj^, /3=/3o ~ l=kzi/kxi, v='Ci/Cs, ^=kxj^Cs/ujp — Xoekxi, 
^=ie(f)/Te, Ni^e=i'ni^Jno. Here a;p=(47rnoe^/mi)^/^ is the ion plasma frequency. Note that 
}C=\k{ti)\/kx^ = (1 +/3^(r) + 7^)^/^ is the modulus of the dimensionless and time-dependent 
wave number vector of SFH. 

Note that the time-dependence of /3(t) and /C(r) is provoked by the temporal variability 
of the SFH wave number component along the flow shear: 

ky{T)=ky^ - Rkx^T = kx^l3{T). (10) 

This process of the "hnear drift" in the k-space [13] has the crucial role for the temporal 
evolution of SFH. 

The total energy of perturbations consists of the kinetic energy of ions, compressional 
energy associated with electrons (ions does not contribute, because their temperature is zero) 
and the energy of electric field. The phenomenological expression for the spectral density of 
the total energy is equal to: 

e^\{\v.? + \vy? + \v.? + [1 + e^/C^ll^r}- (11) 
The first order ODE for £{t), derived from (6-9) is: 



= -R 



VxVy + 2f /3(t)$2 



(12) 
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and in the "no-shear" {R — 0) hmit it is the conserved quantity. 

It is easy to verify that variables appearing in Eq. (6-9) obey the following remarkable 
algebraic relation: 

(1 + -f'^)vy - P{T){va, + -fv^) = RNi + canst. (13) 

By means of this relation, taking a derivative of eq. (6) and re-arranging terms, we can 
get the following exphcit second order ODE for the function ^=A^j//C: 

^(2) + Q2^r)* = constxfir), (14) 

where 

/(t)= - 2R/JC\ (156) 

It should be noted that the analogous, but more simple, kind of inhomogeneous second- 
order ODE was derived and analyzed in [12]. The general solution of (14) is the sum of 
its special solution and the general solution of the corresponding homogeneous {const — 0) 
equation: "if = "i/h + ^^g. When i^(r) depends on r adiabatically, implying [12,17]: 

\n{Tf^\<^n\T), (16) 

then the homogeneous equation can be solved approximately. 

For the flows with R<^1 the condition (16) holds for the wide range of possible values of 
|/9(r)|. In other words, since the temporal variability of \P{t)\ is determined by the "linear 
drift" of SFH, (16) is valid at all stages of the evolution of the SFH. When the condition 
(16) holds, the approximate expression for "ifh may be written in the following way: 

C 

'lfh(r)^^=exp[i(^(T) + (^o)], (17) 



where (/p(t) = jQQ{T')dT' . 

As regards the special solution of inhomogeneous equation (14), it can also be derived 
owed to the smallness of the R parameter. In particular, the solution may be expressed by 



the following series [18]: 



^s{r) = const X J^R^'^y^ir), 



n=0 



yo(r) = /(r)/Q^(r), 

1 d'^Vn-l 



(18a) 
(186) 
(18c) 



Since R<^1, the terms with higher powers of R are negligible and the full approximate 
solution of the inhomogeneous equation (14) may be written explicitly as: 



C r-/ / X M constxfir) 



(19) 



When C/const<^l the SFH may be treated as mainly incompressible and vortical per- 
turbation, while when C / const^l it is mainly of the sound- type. 

Having in hands the solution of (14) we are, certainly, able to find all other variables of 
the problem. They may be calculated by ^ and ^^^^ in the following way: 



A^e = /c*/(i + r^ ), 



(20) 



y ^3 
1 



(21) 
(22) 
(23) 



;i + 72)x:2*w _ + ^2 ^ x:^)^ _ constxpic , 

Note that (20-23) are exact expressions, valid for arbitrary values of the shear parameter 
R. For large enough values of R, when adiabatic solution (19) is no longer vahd, Eqs. (20- 
23) may still be used for reproduction of the variables through ^ and obtained, this 
time, by the direct numerical solution of equation (14). 

Below we shall focus our attention on the behavior of 2D perturbations in the XOY plane 
(7 = 0), when the shear parameter R^l. This case admits simple analytical examination 
and exposes soundly the qualitative novelty of the problem. Using expression (11) for the 
spectral energy density £{t) and (20-23), we get for £{t) the following simple expression: 



^(r)^2 



const 



(24) 



The spatial characteristics of the SFH {kx^, ky{T)) and the value of the shear parameter 
R manage the evolution of the frequency of oscillations and the actual intensity of the energy 
exchange between the SFH and the background flow. In particular, temporal variability of 
these processes are essentially induced by the "linear drift" of SFH in the k-space [13,12]. 

For the sound-type {C^O, const — 0) perturbations, as it is evident from (15a), the 
frequency of oscillations varies with the variation of /C(r). Originally, at moderate values 
of /C(t), due to the smallness of ^, the oscillation mode may be treated as the ion-sound 
wave (f2(r)~/C^(r)). Afterwards, when /C(t) reaches large enough values, the dispersive 
influence of the denominator in the first term of (15a) becomes more and more imperative 
and when ^^/C^(t)»1 the frequency Q(t) already exhibits ion plasma oscillations. Following, 
according to (24), the evolution of energy of this mode (£(T)~n(T)), we find that initially, 
for (3o > 0, at < r < r*=|/3o|/-R, the energy decreases and reaches its minimum at r = r*. 
A while later, it begins to increase at r* < r < oo, when the SFH "emerges" into the 
area of k-space in which ky{T)kxi < (the "growth area" for the sound-type perturbations 
[12]). If the SFH is in the "growth area" from the beginning (/3o < 0), its energy increases 
monotonically. When ^/C(t)>1 the rate of the energy increase becomes less and less and 
the energy asymptotically tends to the constant value. 

When C = and const^O the SFH may be treated as mainly incompressible and 
vortical perturbations. In this case ^c:iconstxf{T)/0,'^{T), while VyC^iconst / JO (t) and 
— constxl3 /10{t). The spectral energy of SFH varies as £^(t)~/C~^(t) and reduces to 
the well-known expression, describing the "transient" growth of the energy of SFH [6,7,13]. 
Transient increase of the energy takes place if initially ky^/k^^ > {(3q > 0) and occurs 
nearby the t^=(5q/R moment of time, when (5{t) tends to zero and /C(t) attains its mini- 
mum value. Such is the behavior of 2D vortical perturbations. One should expect that the 
evolution of three-dimensional (3D) vortical perturbations should be similar to the behavior 
of the analogous structures in incompressible, inviscid fluids, extensively studied in [9] . As 
it was shown in [9], the energy of 3D vortical perturbations grows also noncxponentially, 
but unlike transiently growing 2D perturbations, the energy of 3D perturbations saturates, 
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attaining in asimptotics some constant value. 

Above noted similarity of the behavior of vortical perturbations with the same process in 
usual fluids [12], holds only in the low-R (i?<^l) range. For larger values of the R parameter 
one should expect notable differences between the behavior of vortical perturbations 
in neutral fluids and electron-ion plasma. 

Certainly, in the general case {Cc:iconst), the "vortical" and the "sound- type" evolution 
of perturbations are superimposed on one another. 

Thus, we see that the dispersive nature of plasma, arising due to the violation of the 
quasineutrality for density perturbations, leads to the smooth transition of ion-sound waves, 
amplifled by the process of a shear flow energy extraction, into the ion plasma oscillations. 
As regards to the another type of dispersion properties of ion oscillations, related to the 
finiteness of ion temperature (and, for simplicity, neglected in the above consideration) 
they should lead to the transfer of the energy to the shorter wave length region. The 
characteristic phase velocity of perturbations in the course of the "linear drift" of SFH, 
tends to ion-thermal velocity. Due to the resonant interaction with ions these perturbations 
will be damped (Landau damping) transferring their energy to the particles and accelerating 
them. However, we should notice that such channel of dissipation and coUisionless heating 
of plasma ion component is relevant for perturbations with characteristic time scale smaller, 
than the time scale for electron-ion collisions. The energy can also become one of the main 
sources for the onset of plasma turbulence in such flows. 

Surely, the problem under consideration is quite complex. The comprehensive analysis 
of all possible regimes is beyond the scope of this report and will be published elsewhere. 
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